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Abstract. Two-body energy exchange between stars orbiting massive black 
holes (MBHs) leads to the formation of a power-law density distribution n(r) oc 
r~ a that diverges towards the MBH. For a single mass population, a — 7/4 
and the flow of stars is much less than N(< r)/t r (enclosed number of stars 
per relaxation time). This "zero- flow" solution is maintained for a multi-mass 
system for moderate mass ratios or systems where there are many heavy stars, 
and slopes of3/2<a<2 are reached, with steeper slopes for the more massive 
stars. If the heavy stars are rare and massive however, the zero-flow limit breaks 
down and much steeper distributions are obtained. 

We discuss the physics driving mass-segregation with the use of Fokker- 
Planck calculations, and show that steady state is reached in 0.2 — 0.3 t r . Since 
the relaxation time in the Galactic centre (GC) is t r ~ 2 — 3 x 10 10 yr, a cusp 
should form in less than a Hubble time. The absence of a visible cusp of old 
stars in the GC poses a challenge to these models, suggesting that processes 
other than two-body relaxation have played a role. We discuss astrophysical 
processes within the GC that depend crucially on the details of the stellar cusp. 



1. Introduction 

1.1. Historical overview of the theory 

The formatio n of a stellar cu sp of stars near massive black holes (MBHs) was first 
discussed by iPeeblei ( 19721 ). who argued that a Maxwellian distribution is not 



feasible because stars would be destroyed at unphysical high rates close to the 
MBH. iBahcall fc Woil (jl976l ) performed a more detailed calculation and showed 



that indeed a cusp forms, but with a shallower slope than found by Peebles. The 
Bahcall & Wolf profile of n(r) oc r~ 7 / 4 , where n(r) is the density at radius r, 
has since been confirmed in many theoretical pape r s with different met hods, in- 
cluding Fo kker-Planck (e.g.lCohn fc Kulsrudlll978l: iMurphy et al.l [i99ll). Monte 
Carl o (e.g. Shapiro Marchantlll978l; iFreitag fc Benal2002l ) and iV-body meth- 



ods (jPreto et al.ll2004l : iBaumgardt et al.1 l2004h . and can be understood (with 



hinds ight) by simple dimensional arguments (e.g. lBinnev fe Tremaindl2008l : ISari fc Goldreich 
20061 ) 



In a second paper. [Bahcall &; Wolj ( 19771 ) considered the evolution of stellar 



systems with several masses (denoted with M) and showed that again the re- 
sulting steady state distribution functions fu(E) oc E PM are approximate power 
laws, with the more massive stars having steeper slopes. They found that the 
stellar current into the MBH, Qm(E), is very small ("zero-flow solution") which 
leads to a specific relation between the stellar mass and the logarithmic slope 
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of the DF0, p M = M L /AM H . In the Keplerian limit near the MBH, these DFs 
correspond to power-law density cusps, n(r) oc r~ aM with olm = 3/2 +pm. 
In their analysis, they considered quite moderate number and mass-ratios, and 
concluded that the stars with the lowest masses have slopes oil, > 3/2, whereas 
the stars with the highest masses have slopes an ;$ 2. In contrast to their work 
on single mass-systems, for several decades surprisingly few studies pursued the 
dyn amics of multi-mass systems. 

iMurphv et al.l ( 1 99 ll ) performed a Fokker-Planck study with a mass-spectrum, 
in addition to other physical processes such as stellar collisions, tidal disruptions 
and stellar evolution. More recently, motivated by observations of the stellar 
population surrounding the MBH in Galactic centre (GC), there have been sev- 
eral stu dies of mass-segregatio n near MBHs. T he first applications to the GC 
were bv iFreitag fc Bend (|2002h : iFreitad (120031 1: iFreitae et all (|2006l ). who used 
Henon-typ e Monte Carlo s i mulat ions. These simulations broadly confirmed the 
results by iBahcall Wolf ([ 19771 ) , and showed that cusps should indeed form 
near the MBH in our own GC. 

The planning of the Laser Interferometer Space Antenna (LIS A) formed an- 

other motivation for the study of mass-segregation near MBHs (e.g.. lSigurdsson &; Reed 
19971 1. The reason for this is as follows. When compact stellar remnants spiral 



into MBHs of ~ 10 6 Mq, they emit gravitational waves with frequencies and 
power such that LISA should be able t o obs erve them out to redshifts z > 1. 
It was shown by iHopman Alexander! (|2005l ) that stars can only spiral in suc- 
cessfully if they start very close (< O.Olpc) to the MBH, for otherwise they are 
scattered away from their inspiral orbit. The implication is that inspiral rates 
are very sensitive to the details of the stellar distribution very close to the MBH, 
which is in turn dependent on the amount of mass segregation in the cusp. As 
a result, inspiral rates of stellar black holes (BHs) m ay be higher than those 
of white dwarfs (WDs) (jHopman Alexander! l2006bl ) : even though for typical 
initial mass functions there are many more WDs than BHs in an evolved stellar 
population, mass-segregation drives the BHs to such orbits where they can spiral 
in, such that within ~ O.Olpc there may in fact be more BHs than WDs. In 
addition to inspiral stars, LISA may also ob serve gravitational wave bursts from 
stars during a single fly- by in our own GC ( Rubbo et al. 2006). These stars are 



usually on very short period (~ lyr) orbits, and mass-segregation makes it again 
more likely to observe bursts from BHs than from other stars (jHopman et al. 
120071 1 . 

The interest in mass segregation near MBHs, driven by observ ation, has also 
led to new theoretical developmen ts. lAlexander Hopman! ( 20091 ) extended the 
analysis bv lBahcall Sz Woll (Il977l) to a m uch larger parameter space, and showed 
that the study bv lBahcall Wolj ()1977l ) happened to be at a corner of this space 
where the behaviour is very stable, but that for systems where the massive stars 
are less numerous the evolution is very different; see £|2T1 and Figure dTJ. More 
specifically, they found that the most massive stars ca n have power law slope s 
that are steeper than a = 2, which was ruled out bv IBahcall WoU (|l977l 1. 
This kind of "strong mass segregation" was confirmed in A-body simulations 



1 Throughout this paper the subscripts L, H refer to either the light or the heavy component of 
a two-mass system. 
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bv lPreto Amaro-Seoanel ( 201Cll ). An analyti cal study of mass se gregation with 
continuous mass functions was performed bv lKeshet et al.l (2009), who showed 
that, in the maximally steep limit, heavy stars can develop a cusp as steep as 
n(r) oc r~ 3 . 

In this work we revisit the evolution of a cusp from a duo-mass s tellar 
population which undergoes mass segregation ([Alexander &: Hopmanll2009i ). de- 
termine the time scale for arriving at a steady state distribution, and discuss 
the results in terms of the GC. 



2. Duo-mass systems and strong mass-segregation 



Alexander Hopman (2009) use Fokker-Planck simulations to find the the ap- 
proximate steady state distribution functions (DFs) of a non-evolving spherically- 
distributed duo-mass stellar population around a fixed MBH. Mass is measured 
in units of the mass of a reference star M+, specific energy in units of its velocity 
dispersion e* = cr 2 ., and time in units of its two-body relaxation time at the radius 
of influence, 

3(2vr^) 3 / 2 (1) 



U 



327T 2 'GPMin* In A 



where the Coulomb term is estimated as A = M 9 /M+. Phase space density 
is expressed in units of /* = n*/(27r<7 2 ) 3//2 and distance in units of the MBH 
radius of influence r* = GM m ja 2 . In these units, the dimensionless specific 
orbital energy is defined as x = e/o~* = r*/(2a) (a is the semi-major axis), the 
dimensionless time is defined as r=t/t ir , and the dimensionless DF of each stellar 
mass group as Qm = Jm I '/*■ The evolution of the dimensionless DF, gw, is given 
by the time and energy dependent particle conservation equation that describes 
the two-body diffusion of stars in energy from a fixed unbound reservoir into the 
MBH sink, 

d r/n d 



-gM(x,r) = -x 5/2 —Q m (x,t) . 



(2) 



Or dx 
Qm is the energy fl ow integral which expresses the change in energy due to 



L}m is tne energy now integral wm cn expresses 
two-body scattering (jBahcall h Wolflll979 . fl977h . 
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Equation ([2]) is integrated in time from an arbitrary initial DF until steady state 
is achieved, subject to the boundary conditions that the DF falls to zero at some 
very high energy xrj where the stars are destroyed, and that the unbound stars 
are replenished from a Maxwellian reservoir, 



g M (x > x D ) = 0, g M (x < 0) = C M exg[(al / o- 2 M )x\ 



(4) 



where the constant Cm is related to Sm (the asymptotic number density ra- 
tio of star M relative to the reference star; 5* = 1 by definition) by Cm = 
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(ov/cr^) 3 Sm ■ They assume violent relaxation boundary conditions such that 
Cm = Sm = £*; the steady state DFs do not depend st r ongly on these specific 
choic es of boundary conditions (Ba hcall &: WolJ Il977l ; lAlexander &: Hopman 



20091 ). 

These calculations are used to explore steady state solutions for stellar dis- 
tributions in both the weak and strong mass segregation regime, the nature of 
which they notate using the relaxational coupling parameter A. This parameter 
describes the competition between the self-coupling of the heavy stars and the 
light-heavy coupling in terms of the mass and number ratios, 



N H M 2 H 4 
N L Ml X 3 + Mjj /Ml ' 



(5) 



We reproduce some results of their simulations in Figure ([I]) which shows the 
local logarithmic slopes of the DFs, pl and pn, at x = 10 (which corresponds 
to orbits with a semi major axis of a = O.lpc for r t = 2 pc in the GC) as a 
function of A. For comparison, they take the same mass ratios as modeled by 
BW77, M H /M l = 1.5,3,10. 




Figure 1. The power- law indices Pl,h (thin and thick lines respectively), 
as derived from the logarithmic slopes of <7l,h at x = 10 as a function of A 
for different mass r atios. The logarithmic slopes of the DFs as calculated 
bv iBahcall fc~Wol dl97l are over-plotted as open circles. The transition 
between the weak and strong mass segregation solutions at A ~ 1, which 
is a reflection of the breakdown of the zero-flow assumption as A — > 0, is 
marked. Figure reproduced from I Alexander & Hovmari by permission 

of the AAS. 



Mass Segregation 



5 



In the weak segregation limit (the Bahcall-Wolf solution), A — > oo, which 
is the zero-flow (Qm ~^ 0) limit, the heavy stars dominate the population and 
relax to the single mass cusp an = 7/4 (jpjj = 1/4). The light stars heat by scat- 
tering against the effectively infinite reservoir of heavy stars and diffuse to lower 
energies, thereby settling to a flatter cusp with ai — > 3/2 (pi = M^/AMh — > 0). 
In the strong segregation limit, A — > 0, the light stars, which dominate the pop- 
ulation, behave as a single mass population with oil = 7/4 (pi = 1/4). The 
rare heavy stars sink to the centre by dynamical friction against the effectively 
infinite reservoir of light stars. For low mass ratios, Mh/Ml < 4, where dy- 
namical friction is less efficient, the heavy stars approximately obey the BW77 
relation, p H = (M H /M L )p L = M H /^M L . For higher mass ratios, M H /M L > 4, 
the heavy stars approach the dynamical friction limit, pu — > 5/4. 



The present-day mass function of evolved stellar populations (coeval or 
continuously star forming) with a universal initial mass function, separates into 
two distinct mass scales, ~ 1M© of main sequence stars, dwarfs and neutron 
stars, and ~ IOMq of BHs. In this s ense such systems can be ap proximated 
by a duo- mass system. Furthermore, Alexander Hopman ( 20091 ) show that 
A < 0.1 so that this system should be strongly mass segregated. In the next 
section we calculate the time scale for such a system to form a mass segregated 
steady state cusp. 



3. Time-dependence of cusp formation 

The amount of time it takes to form a stellar cusp around a MBH in a galactic 
nucleus is important for many reasons including, but not limited to, event rates of 
stellar collisions, tidal disruptions of stars and inspirals of compact objects with 
the emission of gravitational waves. In particular, if this time is greater than a 
Hubble time, event rates will be lower than previously predicted. To determine 
this time scale we use Fokker-Planck simulations as described in §211 with initial 
DFs gM{x) oc x -0 ' 5 . We find that steady state is reached in r = 0.2 — 0.3; see 
Figure ©• 

Our result is in agreement with iPreto &: Amaro-Seoand ( 20101 ) who find, 



using both Fokker-Planck and iV-body simulations, that steady state is reached 
in t ~ (0.1 — 0.2) T r i x (r/j), where T T \ x (rh) is the value of the relaxation time 
at the radius of influence of a MBH, and roughly equivalent to t± in Equation 
CO). Values of the relaxation time for MBHs with M, < 6 x 1O 6 M vary in the 
literature, but are generally constrained to within 0.5 — 3 x 10 10 yr. Our results 
imply that quasi-steady, mass segregated stellar cusps should be common around 
MBHs in galactic nuclei fo r this mass reg ime. Hence, assuming GC parameters 
that we observe today (see Merritt 2009), a cusp should form within 



^cusp — 0'2 0.3 t+ 



0.5 - 6 x 10 y yr. (6) 
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Figure 2. Evolution of the distribution function, shown at dimensionless 
times t = {0.1, 0.2, 0.3} with increasing thickness for increasing time. Steady 
state is reached in about r = 0.2 — 0.3. Initially, all DFs were g(x) oc x~ 5 . 
Interestingly, for short times the distribution of white dwarfs (WD), main 
sequence stars (MS) and neutron stars (NS) grows steeper than the final 
distribution which is approximately energy independent. The reason is that 
it takes time for the stellar black holes (BH) to grow a steep cusp and become 
the dominant species. 



4. Applications 



4.1. Galactic centre cusp 



Several studies (e.g. lAlexanderl Il999l : iGenzel et aD 120031 : ISchodel et all 120071 ) 



have shown that indeed a cusp of stars is present in the GC. However, recently 
it has become clear that even though the young s tars do have a power l a w profile, 
there is a dearth of gian t stars within ~ O.lpc ( Buchholz et al. 20091 ; iDo et all 
2009; iBartko et aDl2010l ) . Since the O /B stars are too young to have relaxed by 
the two-body relaxation process discussed here, the conclusion is that currently 
there is no observational support for cusp formation or mass-segregation in the 
GC. The reason for the absence of a cusp is unclear. One possibility is that 
physics not considered here depletes the tig htly bound orbit s of old stars, for 



example as a result of stellar collisions (e.g., Dale et al.l 12009). 



Another option was considered by iMerrittl ( 20091 ). who suggested that re- 



laxational physics does in fact capture the main mechanism for cusp formation, 
but that the age of the GC is too short compared to the time scale for a cusp 
to have formed today. His estimate of the rela x ation time is > 2 — 3 x 10 10 yr, 
and in his Fokker-Planck simulations, IMerrittl ( 20091 ) does not find the forma- 
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tion of a relaxed cusp within that time. This is in contrast to the result we 
found in £|3"! where a c usp f orms in 0.2 — 0.3 t*, similar to the findings by 



Preto fc Amaro-Seoane (120101 ) . Thus, even for a relaxation time of 30 Gyr, the 
central pc would form a cusp in less than a Hubble time, starting from the cur- 
rent situation. The reason for this discrepancy in the time scale for reaching 
a steady state solution is not yet clear, but ma y be related to the back reac - 
tion of the BHs on the light main sequence stars (jPreto &: Amaro-Seoane 2010). 
We note that our findings do not imply that the GC is necessarily relaxed and 
formed a cusp: it is possible that in the past the density at the radius of in- 
fluence was lower, that it has taken a Hubble time to contract^ to the current 
configura tion, and tha t the remaining time to form a cusp is less than a Hub- 
ble time (|Merrittll2009l ). However, this argument appears to be more attractive 
when the remaining time for cusp formation exceeds a Hubble time. 

4.2. S-stars 

Relaxational dynamics cannot account for the presence of a cluster of B-stars 
known as the "S-stars" in the inner ~ 0.03pc near the MBH in the GC (ISchodel et al 
20021 : iGhez et al.ll2005l : lEisenhauer et al.ll2005l : iGillessen et ail 120091 ) . since the 
time scales are too long for young stars (with ages between ~ 20 — 100 Myr) 
to diffuse to such distances. There have been many suggestions of how the 
S-stars have reached their current orbits. Currently the most promising sce- 
narios seem to be those where a binary star is tidally disrupted by the MBH 
(|Hillsl 119881 : lYu Tremainel 120031 1 . Such binaries may originate very far from 
the MBH, where th ey are perhaps dri ven to eccentric orbits by triaxiality or 
massive perturbers ( Perets et al.l 120071 ) . Alternatively, binaries which form in 
eccentric accretion discs can be pumped to very high eccentricities by a gravi- 
tational instability ( Madigan et al. 20091 ). In either CclSG, Si B-star may end up 
on a tightly bound orbit, but the orbit will be very eccentric (e~0.99). Post- 
capture dynamics must then account for the mapping of the initial, very high 
eccentricity distribution to the current distribution in the GC. 

The mechanism res ponsible for changing the eccentricities is most likely 
resonant relaxation (RR; [Rauch fc Tremaind Il996l ) , a secular process that very 
efficiently randomizes the angular momenta of orbits if the precession time is 
much longer than the orbital time-scale. For the S-stars, the resonant relaxation 
time is of order 

T RR ~100Myr^, (7) 

where M* is the typical mass of stars. For M* = IOMq, appropriate for 
mass segregated BHs, this time-scale is s hort enough to randomize the orbits 
( Hopman Alexander! l2006al : iLevhl [2007). Post-capture evolution was studied 
m detail by|Perets| |2009) with the use of A-body simulations. These simulations 
do not include general relativistic precession, which is potentially important for 
the highly eccentric orbits that are considered here, but they clearly show that 



2 Our treatment of the Fokker-Planck equations does not evolve the stars that are unbound 
to the MBH, such that we cannot follow the contraction of the cusp. For a study that 
includes self-consistent mass- segregation and evolution of the ambient cluster of stars, see 
iPreto fc Amaro-Seoand (pOlCD . 
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the distribution evolves on time-scales of the order Trr. For M* = IMq however, 
this time scale is too large for the S-star orbits to have evolved from high ec- 
centricities to their currently observed orbits. Only if there are mass-segregated 
BHs within the S-star radii can RR explain their orbital parameters. 

4.3. Gravitational waves 

As stated in the introduction, most compact remnants that spiral in successfully 
onto MBHs to become detectable gravitational wave sources originate from dis- 
tances very close (< O.Olpc) to the MBH. It is therefore important to what 
extent such orbits are populated, and mass-segregation helps to increase this 
number. In spite of this, even if there is no cusp in the GC, this would not 
necessarily imply that other, similar galaxies do not have cusps either. If, as we 
find, the time for cusp formation from the current state is less than a Hubble 
time, then one would expect that most galaxies with MBHs similar to the GC 
did in fact form cusps, even if our Galaxy di d not. Furthermore, most inspi- 
ral sources originate from lower mass MBHs ( Ho pmanl 120091 : iGairl l2009h : such 



galaxies have even shorter relaxation times, and are thus more likely to have 
formed cusps. The predicted number of detectable LISA sources will therefore 
not depend strongly on the question of whether a cusp is present or not in the 
GC. Co nclusions regarding the absence of a gravitational wave background from 
fly-bys ( Toonen et al J [20091 ) will likewise not be strongly affected. 
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